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Abstract 

We study the dynamics of the striped phase, which has previously been suggested 
to be the ground state of a doped antiferromagnet. Starting from the t — J model, 
we derive the classical equation governing the motion of the charged wall by using 
a ficticious spin model as an intermediate step. A wave-like equation of motion is 
obtained and the wall elasticity and mass density constants are derived in terms of 
the t and J parameters. The wall is then regarded as an elastic string which will be 
trapped by the pinning potential produced by randomly distributed impurities. We 
evaluate the pinning potential and estimate the threshold electric field which has to 
be applied to the system in order to release the walls. Besides, the dynamics of the 
stripe in the presence of a bias field below the threshold is considered and the highl- 
and low-temperature relaxation rates are derived. 
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The discovery of a deep connection between superconductivity and quantum antiferro- 
magnetism in the phase diagram of the cuprate perovskites has stimulated various attempts 
to understand the effects of dilute holes in a spin 1/2 Heisenberg antiferro magnet. The prob- 
lem has usually been addressed by assuming that a doped antiferromagnet can be described 
by a gas of holes with uniform density. However, several calculations ]^J^,[|,£j^3||7J suggest 
that in this system there is a modulation of the charge and spin densities, i.e., the holes 
cluster along lines which separate undoped antiferromagnetic domains (striped phase). 

Experimentally, recent measurements also indicate that striped order indeed occurs in 
doped planar antiferromagnets. In the insulating nickel-oxides, such stripe modulations were 
reported [§,[|,|10| and the data were consistent with multiband Hubbard model calculations 
|TT| . For the case of the copper-oxides, elastic neutron diffraction experiments have revealed 



static stripe order in the non-superconducting compound Lai^gNdcuSro.i^CuC^ [12| while in 
the superconducting compound La2_ x Sr x Cu04 inelastic scattering peaks at incommensurate 
wave vectors suggest the existence of a very similar, albeit slowly fluctuating, striped phase 



Although incommensurate spin fluctuations are not observed in the low doping re- 



gion of the cuprates, muon spin resonance and nuclear quadrupole resonance experiments 
on La 2 _xSr x Cu0 4 with < x < 0.018 [T5|] have been successfully interpreted within models 



that presume a striped structure JT6f. The surprisingly strong suppression of superconduc- 



tivity in these materials for even low Zn doping was also found to be consistent with the 



existence of stripes |L7 



The objective of the present work is to study the dynamics of this striped phase in the low 
doping regime, where interactions between neighbouring stripes is assumed to be negligible. 
We treat the problem on the basis of the t — J model and establish the connection to the 
discrete elastic theory. The domain wall is then considered from a phenomenological point 
of view, i.e., as an elastic line trapped by the pinning potential produced by impurities. The 
depinning of the line from the potential well by applying an electrical field perpendicular to 
the stripe structure is investigated, and the threshold field corresponding to the onset of a 
state with mobile lines is determined. Finally, the relaxation process in the presence of a 



bias field below the threshold is considered and the classical and quantum decay rates from 
the metastable state are computed. 

The first calculations of the striped phase in 2D antiferromagnets with holes have been 
done within the Hubbard model in the vicinity of half filling. They were based on the self- 
consistent Hartree-Fock formalism and have been performed for small and large values of the 
ratio U/t. The results can be summarized as follows: In the small U -approximation, vertical 
domain walls (parallel to the x- or y-axis) are stable whereas for large U diagonal walls 
are energetically more favorable ||||. This crossover from vertical to diagonal stripes was 
numerically calculated to happen at U/t ~ 3.6 (see Fig. 1). 

In the U/t ^> 1 limit, the Hubbard model with almost half-filled band can be reduced 
to the t — J model with effective exchange constant J = 4t 2 /U Jl8|]. Contrary to the 
Hubbard model, where each site of the lattice corresponds to four possible states, in the 
t — J model only three states are allowed, since double occupancy is forbidden. Hence, the 
dimensionality of the Hilbert space in the latter model is much smaller and for finite clusters 
an exact diagonalization of the system, without using the Hartree-Fock approximation, is 
possible. These studies of the t — J model have been recently performed and they confirm 
the results obtained previously: by exactly diagonalizing small systems, Prelovsek and Zotos 
have verified that a striped phase with holes forming domain walls along the (0,1) or (1,0) 
direction arises for J > J s ~ 1.5t, while in the regime 0.4 < J/t < 1.2 the domain walls 
appear along the (1,1) direction. Besides, they found signs of phase separation into a hole- 
rich and a hole-free phase at J > J* ~ 2.51 In Fig. 1 we grafically represent both, the results 
obtained from the t — J and from the Hubbard model in order to make their comparison 
clearer. 

The obtained results for the striped phase can be understood on simple physical grounds. 
In Figs. 2 (a) and (b), diagonal and vertical stripes (ensemble of holes) embbeded in a AF 
background are respectively represented. The holes can jump from one site to another and 
during this process they gain kinetic energy t. An inspection of Fig. 2 leads us to conclude 
that the diagonal configuration favours the dynamics (in this case the holes can move in 
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2D, horizontally or vertically), whereas for the vertical stripe the holes are confined to move 
only in the horizontal direction. On the other hand, in the vertical configuration the holes 
are closer to each other and hence there is a gain in the segregation energy J. Then the 
final behavior results from the competition between the kinetic (t) and the segregation (J) 
energies (see Fig. 2): When t > J = 4t 2 /U, i.e., for U/t large, the kinetic energy dominates 
over the segregation energy and the stripes are diagonally arranged. When the gain in the 
segregation energy becomes more relevant, for t < J, the vertical (horizontal) formation 
arises. 

Recently, the stripe dynamics has been studied from a different perspective: the domain 



walls were regarded from a phenomenological point of view as vibrating strings The 
motion of a single stripe was initially considered, and the analysis was then extended to 
the more complicated regime involving many (interacting) domain walls. Here, we adopt a 
similar description. However, instead of using the initial assumption that the stripe can be 
described by an elastic line, we start from the t— J model and show that the classical equation 
of motion describing the stripe dynamics in the limit of long wavelength displacements has 
a wave-like form. In this way, we deduce the phenomenological mass and elastic coefficients 
for the string in connection with the "microscopic" t and J parameters. For the sake of 
simplicity, we concentrate on the vertical configuration. This implies that our considerations 
hold for J* > J > J s , i.e., 2.5t > J > 1.5t JF|. Besides, we study the problem in the dilute 
limit, when the doping concentration is low and we can investigate the behavior of a single 
stripe (chain of holes). 

In order to derive the quantum equation of motion describing the dynamics of the n-th 
hole in the chain, we map the initial t — J Hamiltonian onto a quantum spin chain problem 
with large spin. Replacing the spin-chain operators in the quantum equation of motion 
by their classical values and considering the long wavelength limit, a classical wave-type 
equation of motion is derived for the stripe. The justification for this classical approach 
rests on the assumption of the existence of zero mode excitations, i.e., we assume that a 
continuum description of the problem is possible and that the discreteness of the original 



lattice formulation is irrelevant in the long wavelength limit. Recently, it was shown p0| that 



the striped phase undergoes a roughening transition, i.e., the flat (gapped) phase present at 
higher values of J ft becomes rough (gapless) as J/t is reduced. Hence, our approach applies 
for the rough phase of the stripe, in which the underlying lattice structure is unimportant 
and the behaviour of the system is governed by a gaussian fixed point. 

Since the equation of motion describing the stripe dynamics has a wave-like form, the 
stripe is from then onwards regarded as an elastic string. The random pinning potential 
due to the presence of impurities is evaluated and the activation energy barrier is calculated 
for the case of a bias electrical field applied perpendicularly to the stripe, both in the 
high temperature phase, where quantum fluctuations can be neglected, and in the quantum 
dominated low temperature phase. The threshold field is estimated and we briefly discuss if 
and how stripe depinning could be experimentally observed. Besides, the dissipation effects 
are considered and their influence in the stripe depinning is discussed. 

The outline of the paper is as follows: in section I, the model is introduced and the 
equation of motion is derived. The pinning potential and the threshold field are calculated 
in section II. In section III we investigate the classical and quantum relaxation regimes at 
fields below the threshold. The discussions and conclusions are presented in section IV. 

I. THE MODEL 

The t — J Hamiltonian describing a chain of holes (stripe) embbeded in an antiferromag- 



netic background is 

HtJ = ~t c t c i° + J 4 - \ n i n j) ■ W 

<ij> <ij> ^ ' 

Here, H t j acts in the truncated Hilbert space consisting of states for which the double 
occupancy of any site is forbidden, i.e., rii — 0,1. 

We apply this Hamiltonian for investigating the dynamics of an infinitely long stripe. The 
linear concentration of holes in the stripe is assumed to be one hole per one lattice constant, 
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so we can neglect charge and spin fluctuations within the stripe. The stripe represents a 
domain boundary (DB) dividing the antiferromagnetic plane into two Neel's phases. We 
consider these two phases to have opposite staggered magnetization (antiphase DB). Hence, 
the string can move without disturbing the initial magnetic order, i.e., it is delocalized (see 
Fig. 3a). This case is opposite to the other hypothetical case, when these two staggered 
fields are equivalent (phase DB) and the string is localized, since its motion is strongly 
frustrated by the surrounding magnetic order (Fig. 3b). 

Here, we concentrate on the dynamics of the vertical configuration and assume that the 
stripe is located along the ^-direction. Further, we make one additional hypothesis, which 
reduces the problem to a directed, but discrete, polymer problem: the holes are constrained 
to move along the transversal ^-direction only, i.e., we neglect overhangs. This assumption 
is true for the shortest excursions of the holes from the initial vertical position of the string 
(up to one step) and therefore seems to be plausible as far as one considers only small 
magnitudes of the oscillations. Hence, the supposition holds for the J > t case, when larger 
excursions of the holes are suppressed by the dominating exchange J— term. In the opposite 
J <t case, the hypothesis ceases to be valid, and the one-dimensional approximation cannot 
be used anymore. This corresponds to the more complicated case of the diagonal DB, which 
will be treated in a separate work. 

Let us denote u n the displacement along the x— direction of the n— th hole from its 
equilibrium position. We describe the relative displacement of two neighbouring holes by 

= u n — u n _i and treat this value as the z— component of some effective local spin a n . 
This spin is related to the n— th segment of the stripe. Displacement of the n— th hole by 
one step along the positive direction of the x— axis increases o z n (a z n — > a* + 1) and decreases 
°n+i(°n+i ~* a n+i ~ !)■ Their sum remains unchanged. Hence, the motion of the n— th hole 
is described by the spin operator — ta~ 2 (a^ +1 a~ + cr~ +1 er+). The a~ 2 factor provides the 
conservation of the norm under the action of this term on the spin-wave function. Besides, 
by accounting for the different Neel order on the left and right hand sides of the stripe, we 
see that any increase of the relative displacement of two neighbouring holes by unity results 
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in a increase of the exchange energy by J/2. Thus, the contribution of the n— th segment 
to the exchange energy of the string is (J/2) | o z n |. 



A similar problem was investigated by Eskes et al. plf , but the Hilbert space in their 
work was restricted in a different way: only configurations where two neighbouring holes are 
separated by one lattice unity were allowed. Hence, they mapped the relative displacement 
coordinate onto a spin one chain and determined the phase diagram features, in the same 



spirit as in Refs. p2| , |23| . Here, we consider a broader Hilbert space and search for the 
classical equations of motion instead of the quantum phase-diagram. 

The projection of the initial t — J Hamiltonian onto the new Hilbert space gives rise to 
the spin-chain problem with the effective Hamiltonian 

H = ~ XX^n+l°n + Vn+lVn ) + i J2 \ a n I ( 2 ) 
° n Z n 

and the standard commutation relations 

The effective spin a is the maximal relative displacement of two neighbouring holes. 

The configuration of the stripe is determined by the z— projection of the spin, a z n . The 
quantum equation governing the motion of the cr* operator is 

h 2 a z n =-[[a z n ,H),H). (3) 

The commutators involving the J— term of the Hamiltonian (which contains the | a z n |) can 
be more easily evaluated with the help of the general relation 



/«J] = oZ{f«) ~ f« ± Km))- (4) 



Here, /(cr^) is an arbitrary function of the operator o z m . 

The next step is to evaluate the commutators on the right-hand side of Eq. (|3|). Then, 
we consider the classical limit o — > oo and replace the spin-operators o z n and by the 
c-numbers: 



— > <rexp(±20 n ) cosa n . 



(5) 



In the long wavelength limit, the discrete variable n can be replaced by a continuous one y. 
Hence, 



Wi - 0n |= a 



dy 



< 1, 



Gin — 0* 



1 du 
ody 



< 1 



(6) 



and the equation of motion (H) acquires the simple wave-like form 

/ d 2 tJa 2 d 2 \ 

[w 2 -irw 2 ) a = - (7) 

Here, t' denotes the time and a is the lattice spacing. A more detailed derivation is given in 
the Appendix A. 

Due to the differential relation between the a and u variables, a = (a/a)(du/dy), the 
displacement of stripe obeys the same wave-like equation, 

/ 8 2 Ua 2 d 2 ^ 



u = 0. 



(8) 



\dt' 2 h 2 dy 2 

This equation describes the long wavelength oscillations of the stripe around the equilibrium 
position u = const. The corresponding action reads 



/OO /"OO 
dy / dt 1 
-oo J 



h 2 (du 



J ( du^ 



(9) 



2ta 3 \dt> J 2a\dy) 

The problem then has been mapped onto a massive string with linear mass density p = h 2 /ta 3 
and elastic tension coefficient C = J/ a. 

Finally, one can observe from Eq. (H) that the long wavelength elementary excitations 
of the string are gapless and have the phonon-like dispersion relation 



uj = ck, 



(10) 



with phase velocity c = aytJ/fi. 

We emphasize again that our derivation assumes the validity of a continuum description. 
Hence it holds for the rough phase. 
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II. THRESHOLD FIELD 



Up to now, we have only considered the dynamics of a stripe embedded in an antifer- 
romagnetic background. However, in real systems, we have also to take into account the 
presence of randomly distributed point-like impurities which act to pin the lines, leading to 
a glassy phase with trapped stripes. In order to reduce the pinning barrier, we apply an 
external electrical field perpendicularly to the stripes formation. Then, we determine the 
threshold field above which the potential barrier vanishes and the stripes can flow through 
the sample, as in a liquid state. 

The free energy describing an elastic string along the y-direction, which tends to move 
due to the action of an externally applied electrical field E competing against the pinning 



barrier V pin is 



/oo 
dy 
-oo 



C (du\ 2 eEu 
TT I t; I +V vin - 



(11) 



2 \dyj ' pm a 

The next step now is the evaluation of the pinning potential V P i n . Let us consider that 
the pinning mechanism is produced by the ionized acceptors which sit on a plane parallel to 
and close to the CuO-plane. An impurity with two dimensional coordinates R produces at 
the position r in the CuO-plane the Coulomb potential G(R — r) 

G(R-r) = -f- (12) 
e\R — r\ 

The total potential felt at position f in the CuO planes can then be written as 

V pin (?) = -J d 2 RAf(R)G(R - f), (13) 

(X J 

where M{R) is the number of impurities (zero or one) at the position R. Hence, M{R) = 
(Af(R)^j . Let us denote the average number of impurities per one site (j\f(R)^ = v, where 
(. . .) represents average over the disordered impurity ensemble and < v < 1. If the 
impurity distribution would be uncorrelated, we could write the density-density correlator 

as 
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Po = (M(R)M(R f )) - (Af(R)) (M(R')) = v(l - u)S(R - Rl). (14) 

— * —* 

Its Fourier-transform is a constant for any value of the wave vector k, po(k) — u(l — u) — 
const. This correlator describes long-wave fluctuations with any k, including the small ones. 
Although the small-/c fluctuations are allowed by the statistics of the completely disordered 
state, in reality they are strongly suppressed by the long-range Coulomb interaction of the 
impurities, since such fluctuations provide large value of the Coulomb energy ~ k^ 1 . Hence, 
the real correlator should be suppressed for small k at some scale k < A. This cut-off can 
be introduced by the following choice of the correlation function: 

v(l — v) , 

= TTaW' (15) 

Since the only length scale in the system of the impurities is the average distance between 
two neighbor impurities u = az/" 1 / 2 , the A parameter should be estimated as A ~ 1/u. 
Further, the two-point correlator of the impurity potential is given by 

= \ I d 2 R d 2 R'G{R')G{R - r)p{R, R'). 
The Fourier transform of /C(f) reads 

/„2\ 2 A „2 



4tt 2 

eaj (PTA 2 )a 2 



K{k) = G\k)p{k) = {v-v>)\-) fu2 ■ (16) 



where G(k) = 2ne 2 / (ek) is the 2D Fourier-transform of the Coulomb potential G(r). Rewrit- 
ing then the correlator of the potential in real space we find 

/C(f) = 2vr (^) (u-u 2 )K (\r), (17) 

where Ko(Xr) is the modified Bessel function. 

Let us now define e P i n = J dyV P i n . If the stripe interacting with the random pinning 
potential is stiff, the average pinning energy (e pin (L)) d of a segment of length L is zero. The 
fluctuations of the pinning energy, however, remain finite, 
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E 2 . 
P m /d 



with 



4 f L 'dydy'1C(0,y-y')~7L (18) 

CI j 



7= a^. (19) 



1/2 

Here, e c = e 2 /ea denotes the Coulomb energy scale. The sublinear growth of (e^^L)^ is 
due to the competition between individual pinning centers. The dynamic approach to this 
problem was introduced by Larkin and Ovchinnikov |3B| in the "Collective Pinning Theory" 
(CPT) for describing the dynamics of weakly pinned vortex-lines in the high temperature 
superconductors. A scaling approach was also considered in connection with the pinning 



problem in charge-density- wave systems [^j,^ . The results of the CPT can be summarized 



as follows ||29|| : Eq. (|18D implies that a stiff stripe is never pinned, since the pinning force 
grows only sublinearly, whereas the electrical driving force increases linearly with length. On 
the other hand, due to the elasticity, the stripe can acommodate to the potential on some 
"collective pinning length" L c . Hence, each segment L c of the stripe is pinned independently 
and the driving force is balanced. 

Our task now is to determine this length L c . The evaluation of the free energy by using 
dimensional estimates provides 

^[«,L]~c£-^-— ■ (20) 



By minimizing !F[u, L]/L with respect to L at zero bias field [2]| we determine the collective 
pinning length L c along the string, 

«7* . 0i Lc j^y\ (21) 



SL 



Assuming that u ~ A 1 ~ a/y/u, the average impurity spacing, and using eq. fllTf ) we find 

L c ^au~^( £ -i) 2/ \ (22) 



where the elastic energy E\ = Ca. Experimentally, it is difficult to measure the collective 
length L c . However, the threshold electric field E c corresponding to the vanishing of the 

11 



barrier is a quantity which can be easily experimentally determined. The critical value E c 



can be estimated by equating the pinning energy y/^yT7 c to the electric energy eEuL c / a. We 
then obtain 

Ec= ±^„^(±y\ (23) 

eu y L c ea KeiJ 

For typical high-T c materials, such as La 2 _j y Sr i/ Cu0 4 , J ~ 0.1 eV, a pa 4 A and e = e ~ 
30 |0|. Hence, one can estimate the elastic energy si = J pa 0.1 eV and the Coulomb energy 
e c = e 2 / (ea) pa 0.1 eV. For a doping concentration in the antiferromagnetic insulating phase, 
for instance, for v = 10~ 3 , we then obtain the collective pinning length L c w 10 3 A and the 
critical electrical field E c pa 10 3 V/cm. 



III. CLASSICAL AND QUANTUM RELAXATION PROCESS 

We have estimated the threshold field E c corresponding to the onset of the motion of 
depinned stripes. Next, we are interested in studying the relaxation process taking place at 
applied fields E < E c . In this case there is a finite pinning barrier preventing the motion 
of the stripe, but it can still jump over (under) the barrier due to thermal (quantum) 
fluctuations. 



A. Classical limit 

At high temperatures T, the decay rate T t is given by the Arrhenius law, T t ~ 
exp(—U/T). The activation energy U can be determined by extremizing the free energy. 
In order words, within the semiclassical approximation, the energy barrier U is nothing but 
the free energy T evaluated at the saddle point configuration u s , U = J-[u s ]. As far as we 
are interested in evaluating the decay rate only within exponential accuracy, we can safely 
neglect the dynamical terms (as for instance, the kinetic one) in the free energy, since these 
terms would give a correction only to the prefactor multiplying the exponential function. 
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By substituting the collective length L c as given by eq. (^) into the free energy (p0[) , we 
can estimate the collective pinning energy barrier 



U c ~ ^L c ~ ^- 1 /6( £ ^2 ) l/3_ (24) 

For the case of L^-t/Sr^CuC^ with z/ = 10~ 3 considered in the previous section, we estimate 
the pinning barrier to be of the order of 10 3 K. The barrier exhibits a weak dependence 
on the doping parameter, U c oc v^ 1 ^ . Notice that this barrier height is compatible with 
the estimates presented in jn| for the binding energy of holes in a domain wall. Besides, 
this value is also comparable to the one obtained for the vortex creep process in high-T c 
superconductors, when U c ~ 10 2 — 10 3 K. 

The dynamics of the stripe in the presence of the pinning potential is simply an example 
of the more general problem of elastic manifolds in quenched random media. Investigations 
of the statistical mechanics of this object have shown that a stripe confined to move in a 
plane is always in a pinned phase, 

u(L)-u(0)] 2 )) ~ ul{^]\ L>L C , (25) 



with a wandering exponent ( = 2/3 |5T , 3"2" |. Here, ((•••)) denotes the full statistical average 



over dynamical variables (thermal) and over disorder, L is the distance along the stripe, 
and u c and L c are transverse and longitudinal scaling parameters, respectively. In our case, 
u c ~ A -1 (scale of the disorder potential) and L c is the collective pinning length. 

Besides, it was also found that competing metastable states that differ from one another 
on a length L are separated by a distance 

/ T \ C 

u{L)~u c [—) , L>L C (26) 



and a typical energy barrier 



U{L)^U C [^] K \ L>L C , (27) 



where U c denotes the scaling parameter for energy. For the single stripe problem, U c reduces 
to the collective pinning energy. The free-energy functional at low driving fields E E c is 
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(28) 



The problem now has been reduced to a nucleation process |33j . If a nucleus with length 
L larger than some optimal length L opt is formed, the system will move to the next minimum. 
On the other hand, if the activated segment is smaller than the optimal one, the nucleus 
will collpase to zero. The optimal nucleus can be found by extremizing the free-energy, 
dLJ~(L)\L=L ovt = and we obtain 

> EA 1/(2-0 



Lopt(E) ~ L c 



E 



(29) 



Inserting back into the free energy ( p5|) we verify that the minimal barrier for creep 
increases algebraically for decreasing bias field, 



' ft \ M 

U(E) ~U C {^ 



(30) 



with fi = (2( — l)/(2 — C) = 1/4. Hence, the system is in a glassy phase, with a diverging 
barrier in the limit of vanishingly small applied electrical fields. 

Another interesting limit to study the dynamical behavior of the stripe, is at fields below 
but close to the critical field E c , i.e., at E c — E <C E c . In this case, the effective potential 



given by the pinning and the bias electrical field terms can be written as [2{J 



V eff(u) 



u v 

UfJ 



u y 

UfJ 



(31) 



with V F ~ V c (l - E/E c fl 2 and u f ~ u c (l — E/Ec) 1 / 2 . The critical potential barrier 
V c = eE c u c /a. The energy of a distortion «j? of the stripe on a scale Lp is estimated to be 



S(u F ,L F ) 



eff 



(u F ) 



(32) 



2 \L F J 

The competition between the barrier to be overcome Vp with the elastic energy density 



Cup/L 2 F determines the length of the saddle-point configuration 



jC ( E\-^ 

Lfs ~ Uc\ 77 77~ 



V, 



E r 



(33) 
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Finally, the energy barrier for thermal activation of the stripe out of the pinning potential 
reads 

, p x 1/2 / 77. s 5/4 / 77. s 5/4 

U(E)~« C V C ( V ) (l-_) ~ 17.(1-^) ■ (34) 

Eqs. (|30|) and (p4[), together with ( ^4|) are the main results of this section. From an ex- 
perimental point of view, it would be easier to observe the creep of the stripe near criticality, 
where the thermal process is described by Eq. ([34]) and the activation barrier ( ^4[) can be 
reduced by one or two orders of magnitude due to the presence of the bias electrical field. 



B. Quantum limit 

At low temperatures, we expect the decay process to be driven by quantum fluctuations. 
In this case, the dynamical terms become essential since they are related to the so-called 
traversal time. The tunneling rate is then T q ~ exp(— B/h), where B is given by the 
Euclidean action Se of the system (the action in an imaginary time formalism) evaluated 
at the saddle point solution, B = Se[u s ]. The total Euclidean action describing the elastic 
domain wall in the presence of random impurities and an external electric field is 



/oo /*oo 
dy / dr 
-00 JO 



p ( du\ 2 C ( du\ 2 Tr eEu 



2 \drj 2 \dyj P " L a 
The tunneling time r c can be estimated by equating the kinetic and elastic terms, 



(35) 



p-f ~ C-f. (36) 



We then obtain r c ~ L C ^J p/C. Substituting r c into the Euclidean action ([35]) , we find 

ST ~ r c L c ^$ ~ ulJ^C ~ - S ■ (37) 



It is important to notice that the extremal value of the action B™ does not depend on the 
collective pinning length L c and hence it is independent of the pinning potential. 

Next, we account for dissipation effects in order to generalize our model. Since we 
cannot determine the friction coefficient 77 from microscopic calculations (we are considering 
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frozen Neel phases for describing the antiferromagnetic background), we will evaluate it 
in a phenomenological way. Moreover, we restrict ourselves to the simplest case of ohmic 



dissipation and study the problem within the framework of the Caldeira-Leggett model ]31 

In the overdamped limit, when ^r c /p <C 1, we can neglect the massive term in equation 
(B5I) and substitute 



P_(du\ ^ roc rn_ 

2\dr Jo 4tt 



u r — u[t 



T — V 



(38) 



The tunneling time now can be obtained by comparing the dissipative and elastic terms. It 
reads ~ rjiL^/C and the corresponding minimal action for the tunneling process is 

B d c ~ r?L c ^f- ~ Vl ulL c . (39) 

Now, we estimate the calculated values in order to better understand our results. The 
friction coefficient per unit length r\\ = rj/a can be estimated from the known data for the 
metallic phase (y = 0.1). By using the Drude formula, we can evalute rj = ne 2 /a and 
the corresponding relaxation time r = m/i]. Here, n is the hole concentration per unitary 
volume Vo, n = u/Vq ~ 10 _1 /200A 3 = 0.5 x 10 21 cm -3 , a is the normal state conductivity, 
a ~ 10 3 1 cnr 1 |^| and m is the effective tight-binding mass of the carriers in the CuO 
plane, m/m e ~ 1 |36|, with m e denoting the free electron mass. We then obtain rj ~ 10~ 13 



g/s and r ~ 10~ 14 s. The relaxation time is assumed to be independent of the doping, and 
to have the samer order of magnitude for the metallic and for the insulating states. The 
tunneling time r c ~ 10 -12 s and t^J ~ 10~ 10 s. This indicates that the quantum dynamics is 
overdamped, and that a proper theory for describing the relaxation process should account 
in a more accurate way for the dissipative term. 

Let us now calculate the correction in the minimal action due to the presence of the 
electrical field. In general, a quantum problem can be regarded as a (d + l)-dimensional 
generalization of a classical ci-dimensional problem, with the imaginary time being considered 
as an additional dimension. However, one should notice that the disorder potential fluctuates 
in space, but not in time. At high temperatures, the stripe jumps over the barrier and the 
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time needed for the jump is irrelevant. At low temperatures, instead, this time is essential. 
Hence, for the classical motion the stripe can choose optimal barriers [see Eq. whereas 
for tunneling the relevant barrier scales like the average barrier U C (L/L C ) (see |29|] for a more 
detailed discussion). As a consequence, we obtain that in the limit of low driving fields 
E -C E c , the minimal action reads 

B - B \-i) > (4°) 

were fi q = (1 + 0/(2 — C) an d B c is B™ or B^ [given by Eqs. (|37|) or (j39|), respectively], 
depending if we are considering the massive or the dissipative limits. Remembering that in 
our case £ = 2/3, we finally obtain a quantum glassy exponent /i g = 5/4. Near criticality, 
the results change considerably: as in the thermal case, in this limit the quantum action 
exhibits a power law behavior, 

fl~«.(l-f)", (41) 
with a = 1 for the massive stripe and a = 3/4 for the overdamped one. 

IV. DISCUSSIONS 

In this work we have succeeded in deducing a wave like equation for the motion of a 
line of holes embedded in a antiferromagnetic phase, which we have chosen to describe by a 
t — J model. This was done through the application of semi-classical methods to a fictious 
spin chain whose local spins can appropriately be related to the displacement difference 
between neighbouring holes along a direction perpendicular to the line itself. Therefore we 
were able to establish a connection between the microscopic parameters of the t — J model 
and the phenomenological mass and elastic coefficients of the continuum theory by using 
the ficticious spin model as an intermediate step. 

We have also extended the well-established "collective pinning theory", which has suc- 
cessfully been applied to vortices in superconductors and charge density waves in quasi 1-D 
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electronic systems, to the case of the wall of holes. In so doing, we were able to estimate 
the threshold field for the depinning of the wall as well as the energy barrier per unit length 
felt by the trapped line. These quantities allowed us to compute the thermal and quantal 
rates for the depinning of a single wall of holes. 

These results can be tested by measuring the conductance of a sample, which presents 
the striped phase in its antiferromagnetic regime, as a function of the external field (voltage). 
From our findings, we expect to have a vanishingly small conductance up to the threshold 
field and then a gradual tendency to recover the ohmic regime for E > E c . This highly 
non-linear behavior of the conductance is analogous to that observed in charge density wave 
systems ||37|| . It is also important to notice that in order to compare our results with the 
experimental ones, one should take into account that inhomogeneities in the barrier height 
distribution may affect the observed threshold field. Actually, the measured value is a lower 
bound for the estimated voltage, and results differing by even one order of magnitude from 
our findings would not be surprising. 

Finally, we would like to say some words about our approach to the damping of the 
wall motion. We have chosen to use an entirely phenomenological approach to the problem 
because it is not possible to describe dissipation from the present microscopic model. The 
first step in this direction would be to allow the spin system to respond to any change of 
the line configuration, i.e., we should assume a finite stiffness for the magnetic system. This 
procedure is planned to be done in a future publication. 

We are indebt with M. Balina, D. Baeriswyl, A. H. Castro Neto, and H. Schmidt for 
fruitfull discussions. This work has been supported by the DAAD-CAPES PROBRAL 
project number 415. N. H. acknowledges financial support from the Gottlieb Daimler- 
und Karl Benz- Stiftung and the Graduiertenkolleg "Physik nanostrukturierter Festkorper" , 
Universitat Hamburg. Y. D. acknowledges financial support from the Otto Benecke-Stiftung. 
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APPENDIX A 

Our aim here is to calculate the commutators on the right-hand side of Eq. With 
the help of the relation (|4]), one obtains 

h 2 a z n =-[[a* n ,H],H] = 

= ^2 | i^n^n+l) (l ^n+l I + I K I ~ I °"n+l ~ 1 I ~ I °"n + 1 l) 
- ( a n °"n+l) (l <+l I + I < I - I <+l + 1 I - K - 1 
+ (°~n°"n-l) (l <-l I + I < I - I <-! - 1 I - I < + 1 



2t_ 2 
V 4 



H r< erf .4.1 (c^ov, + o"„ ^ + cr+cr„ , 5 + a, n a2 



+° r n-l °n + °n + ^n ff n-2 + °n ^n-z) } 

(°"n+l°"n+l + °n+l°"n+l + °"n-l°n-l + ^n-l^n-l + 2o "n-l°n+l + 2o "n-l°"n+l) | 

Now, we consider the problem in the classical limit, when the operators cr* and can 
be replaced by c-numbers. In spherical coordinates, we can write 

a n = a sin ot n , &n = c"exp(±i0 n ) cosa n , (Al) 

where a n = tt/2 — 6 n . Using that 4> n — <p n -i ~ and a n <C 1, we obtain the classical equation 
of motion 

(8t \ 
J +^J K+l ~ 2«„ + On-i) . (A2) 

Further, we take the limit a — > oo and go to the continuous approximation. In this limit, 
we can replace 

a 2 

ot n -> a(y), a n± i -> a(y) ± aa'(y) + ya"(?/), (A3) 

and we eventually obtain the wave-like equation 

/ d 2 tJa 2 d 2 ' 
[dt^ ~ n 2 dy 2 



a = 0. (A4) 
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FIGURE CAPTIONS 

Fig. 1: Regions corresponding to the diagonal and vertical striped phase obtained from 
(a) exact diagonalization in small systems 0; (b) Hartree-Fock calculations ||. 

Fig. 2: (a) Diagonal and (b) vertical stripe configurations. 

Fig. 3: (a) Delocalized stripe (antiphase domain boundary); (b) Localized stripe (hy- 
pothetical case corresponding to a phase domain boundary). 
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